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Abstract 

In this work, the thermal expansion properties of carbon nanotube (CNT)-reinforced nanocomposites with CNT 
content ranging from 1 to 15 wt% were evaluated using a multi-scale numerical approach, in which the effects of 
two parameters, i.e., temperature and CNT content, were investigated extensively. For all CNT contents, the 
obtained results clearly revealed that within a wide low-temperature range (30°C ~ 62°C), thermal contraction is 
observed, while thermal expansion occurs in a high-temperature range (62°C ~ 120°C). It was found that at any 
specified CNT content, the thermal expansion properties vary with temperature - as temperature increases, the 
thermal expansion rate increases linearly. However, at a specified temperature, the absolute value of the thermal 
expansion rate decreases nonlinearly as the CNT content increases. Moreover, the results provided by the present 
multi-scale numerical model were in good agreement with those obtained from the corresponding theoretical 
analyses and experimental measurements in this work, which indicates that this multi-scale numerical approach 
provides a powerful tool to evaluate the thermal expansion properties of any type of CNT/polymer nanocomposites 
and therefore promotes the understanding on the thermal behaviors of CNT/polymer nanocomposites for their 
applications in temperature sensors, nanoelectronics devices, etc. 
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Background 

As technology and modern industry has developed, rein- 
forced composite materials, such as particle- or short- 
fiber-reinforced composites and long-fiber-reinforced or 
sandwich laminates, have been widely applied in the aero- 
space, construction, transportation, machinery, chemical, 
and other industries. In recent years, as a representative of 
new engineering materials, carbon nanotube (CNT) at 
nanoscale has shown superior mechanical, electrical, and 
thermal properties, as well as low density and high aspect 
ratio, which make it an ideal choice for composite 
reinforcement. CNT-reinforced nanocomposite is a multi- 
phase material, and its external macro-physical properties 
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strongly depend on the properties of its constituents and 
complex internal microstructure. Experimental evaluation 
requires large amounts of material samples and a large 
testing work load, giving simulation of the physical proper- 
ties of nanocomposites important engineering significance. 

There has been extensive research on the mechanical, 
thermal, and electrical properties of CNT-reinforced 
nanocomposites. For instance, the thermal properties 
[1-3] and electrical properties of CNT-reinforced nano- 
composites [4,5] have been explored experimentally in 
some previous studies. Moreover, due to the complexity 
and variations of the CNT-reinforced composite micro- 
structure, theoretical analyses and numerical simulation 
methods are common strategies to estimate composite 
physical properties. For instance, diffusion and thermal ex- 
pansion coefficients of CNT-reinforced nanocomposites 
have been studied through micromechanics models with- 
out sufficient atomic scale information [6] or molecular 



Spri 



© 2013 Alamusi et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons 
1"1 Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction 

?S in any medium, provided the original work is properly cited. 



Alamusi et al. Nanoscale Research Letters 2013, 8:15 
http://www.nanoscalereslett.eom/content/8/1 /1 5 



Page 2 of 8 



dynamics (MD) models with very high computational cost 
and complexity [7]. 

In recent years, to deal with the remarkable scale differ- 
ence in CNT-reinforced nanocomposites, multi-scale mod- 
eling has been widely used for predicting the mechanical 
properties [8], electrical properties [9], and thermal con- 
ductivity [10] of the CNT- reinforced nanocomposites. 
However, to the best knowledge of the present authors, 
there has been no report on the multi-scale modeling of 
thermal expansion properties of the CNT-reinforced 
nanocomposites to date. In this work, the thermal expan- 
sion properties of the CNT-reinforced nanocomposites, 
i.e., CNT/epoxy, were evaluated using a sequential multi- 
scale numerical model. The present study focused on the 
effects of two key parameters, i.e., temperature and CNT 
content, on the thermal expansion properties. Moreover, it 
was found that the results of the present multi-scale nu- 
merical model agree very well with those based on theor- 
etical predictions and experimental measurements carried 
out in this work. 

Methods 

To investigate the thermal expansion properties of CNT- 
reinforced nanocomposites, numerical simulations based 
on a sequential multi-scale approach were conducted on 
two types of microstructural models, a uni-directional 
model in which CNTs were uni-directionally aligned 
within epoxy and a multi-directional model in which the 
CNTs were randomly oriented within the epoxy. For 
CNT-reinforced nanocomposites, uni-directionally aligned 
CNTs in matrix can be realized by applying electric [11] 
or magnetic fields during the curing process. The uni- 
directional model was constructed as a two-dimensional 
(2D) axisymmetric model (see Figure 1), and the multi- 
directional model was built up as a 2D plane strain unit 
cell model (see Figure 2). Note that to reduce the compu- 
tational cost, an equivalence conversion principle [12,13] 
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Figure 1 Schematic of uni-directional numerical model, (a) A 

cylindrical model (RVE). (b) Schematic of a quarter axisymmetric 
model. 




from three-dimensional (3D) modeling to 2D modeling 
for short-fiber-reinforced composites was used as a sup- 
porting evidence for the present 2D plane strain multi- 
directional model. 

To construct the sequential multi-scale numerical 
model, we firstly used the axial thermal expansion proper- 
ties of multi-walled carbon nanotube (MWCNT), which 
were obtained from extensive MD simulations at atomic 
scale in the authors' previous work [14]. Secondly, con- 
tinuum mechanics-based microstructural models, i.e., the 
uni-directional and multi-directional ones, were built up 
based on the MWCNT's thermal expansion properties at 
atomic scale and the thermal expansion properties of 
epoxy obtained from experimental thermomechanical 
analysis (TMA) measurements in this work. The detailed 
description of experiments will be provided later. The 
thermal expansion rates s of the present MWCNT and 
epoxy from 30°C to 120°C are shown in Figure 3. As 
shown in [14], the axial thermal expansion rate of 
MWCNT is dominated by MWCNT's inner walls. We 
modeled MWCNT's six innermost walls [14] to obtain the 
approximate axial thermal expansion rate of the present 
MWCNT in Figure 3. 

In the uni-directional and multi- directional models 
used for the finite element analysis, the present multi- 
scale numerical simulations were conducted under the 
following conditions: 

1. The CNT content of CNT/epoxy nanocomposites 
ranged from 1 to 15 wt%. 

2. The length and diameters of the outmost and innermost 
walls of CNT were set as 5 [am, 50 nm, and 5.4 nm, 
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Figure 3 Thermal expansion rates of CNT and epoxy. 



respectively, which are in accordance with the experi- 
mental measurement using a transmission electron 
microscope [9,15]. The properties of MWCNT used in 
the present experiments are shown in Table 1. 

3. We only considered the axial thermal expansion/con- 
traction of MWCNT, and the radial thermal expan- 
sion/contraction was neglected since they are very 
small as identified in [14]. Therefore, CNT thermal 
expansion properties were orthotropic. Other proper- 
ties of CNT were assumed to be isotropic, as well as 
those of epoxy. The detailed material properties in 
simulations are listed in Table 2. 

4. For the uni-directional model, simulations were con- 
ducted using a quarter of the cross section of a cylin- 
der representative volume element (RVE) containing 
a CNT, i.e., an axisymmetrical model (see Figure 1). 
Under thermal loading, some forces along the radial 
direction were imposed on the nodes of the outmost 
lateral surface of the RVE and adjusted through an it- 
erative procedure so that all points on the outmost 
lateral surface moved at the same distance in the ra- 
dial direction to simulate the periodic conditions 
[16]. The length of the polymer was two times longer 
than that of the CNT in Figure 1, implying that the 
short CNTs are distributed evenly in both longitu- 
dinal and lateral directions in a matrix so that the 
RVE is the same for any CNT [16]. 



Table 1 Properties of MWCNT 

Property 



Value 



Fiber diameter (nm) 
Aspect ratio (-) 
Purity (%) 



Average 50 

>100 

>99.5 



Table 2 Material properties 


Property 


CNT 


Epoxy 


Density (g/cm 3 ) 


2.1 


1.1 


Young's modulus (GPa) 


1,000 


3.2 


Poisson's ratio 


0.1 


0.34 


Specific heat (mJ/g-K) 


650 


1,000 


Thermal conductivity (W/mnrvK) 


6.7 


2 x 10" 4 


CTE (K" 1 ) 


From Figure 3 


From Figure 3 



5. For the multi-directional model, there were randomly 
distributed 100 CNTs per model (see Figure 2). This 
model was built up under plane-strain conditions. 
The boundary conditions were applied at the two ex- 
ternal edges which is similar to those for the uni- 
directional model above. In order to reflect the 3D 
characteristics of real nanocomposites, the volume 
fraction should be converted to the half of the real 
one [12,13]. Note that the number of the CNTs in 
this model, i.e., 100, was determined by some trial 
computations, such as testing of models containing 
10, 25, and 50 CNTs. It was found that 100 is the 
minimum number, which can yield isotropic, conver- 
gent, and stable results. This number is just the same 
with that of holes for modeling the effective mechan- 
ical properties of a porous plate [17]. 

Results and discussion 

Uni-directional models 

Firstly, we investigated the influences of temperature and 
CNT content on the thermal expansion properties of 
CNT/epoxy nanocomposites by varying the temperature 
from 30°C to 120°C and CNT content from 1 to 5 wt%. 
The thermal expansion properties vary with temperature, 
as shown in Figure 4. In this figure, the thermal expansion 
rate increases linearly as temperature increases for any 
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Figure 4 Thermal expansion rate of uni-directional CNT/epoxy 
nanocomposite by numerical simulation. 
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loading of CNT. The temperature of zero thermal expan- 
sion rate (or no thermal expansion/contraction) of the 
CNT/epoxy nanocomposites is approximately 62°C, which 
is independent of CNT loading. Moreover, at a specified 
temperature, the absolute value of thermal expansion rate 
decreases with increasing content of CNT. The influence 
of the nonlinear thermal expansion rate of CNT (Figure 3) 
on that of the nanocomposites seems to be small due to 
very low CNT contents in Figure 4. 

Although it is still a technical challenge to uniformly dis- 
perse CNTs for high loading, e.g., over 10 wt%, to numer- 
ically explore the thermal expansion properties in detail, 
the content of CNT was varied from 1 to 15 wt%, and the 
corresponding results are shown in Figure 5 with some 
artificial adjustments due to the big differences in various 
curves. From Figure 5, the thermal expansion rates vary 
nonlinearly with the content of CNT. In the range of 1 to 
5 wt%, the change of thermal expansion rate is obvious. 
Beyond 5 wt%, the increase of CNT content within the 
temperature range (30°C ~ 120°C) results in the absolute 
values of the thermal expansion rate \e\ becoming grad- 
ually smaller and finally converging to a stable value when 
the CNT content reaches 10 wt%. Note that the thermal 
expansion rate is negative at 30°C. 

Multi-directional models 

The ranges of temperature and CNT content in this case 
are identical to those mentioned above for the uni- 
directional models. The variation of thermal expansion 
properties of CNT/epoxy nanocomposites is shown in 
Figure 6 (CNT content from 1 to 5 wt%), in which the 
similar effects of temperature and CNT content are 
observed. In this figure, the thermal expansion rates in- 
crease linearly as the temperature increases for all CNT 
contents. The temperature at zero thermal expansion 
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Figure 6 Thermal expansion rate of multi-directional CNT/ 
epoxy nanocomposite by numerical simulation. 
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rate (or no thermal expansion/contraction) of the CNT/ 
epoxy nanocomposites is approximately 62°C at any 
CNT loading, which is similar to that for the uni- 
directional model. With increasing content of CNT, the 
absolute value of thermal expansion rate decreases. 
Moreover, compared to the uni-directional nanocompo- 
sites (Figure 4), at high temperature, the difference in 
thermal expansion between low CNT content (1 wt%) 
and high CNT content (5 wt%) is much smaller in the 
multi-directional nanocomposites. 

By varying the CNT content from 1 to 15 wt%, the 
obtained results are shown in Figure 7. In this figure, the 
thermal expansion rates vary nonlinearly with the CNT 
content. In the content range of 1 to 5 wt%, the change in 
thermal expansion rate is obvious. Beyond 5 wt% CNT, as 
the CNT content increases, the absolute value of the ther- 
mal expansion rate \s\ becomes smaller gradually. However, 




Loading (wt%) 



Figure 5 Relationship between CNT content and absolute 
value of thermal expansion rate of uni-directional CNT/epoxy 
nanocomposite. (Data of 30°C = Original data x (-2.5); data of 
75°C = Original data x 8). 
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Figure 7 Relationship between CNT content and absolute value 
of thermal expansion rate of multi-directional CNT/epoxy 
nanocomposite. (Data of 30°C = Original data x (-2.5); data of 
75°C = Original data x 8). 
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unlike the uni-directional nanocomposites (Figure 5), the 
thermal expansion rate of the multi-directional nanocom- 
posites still decreases proportionally to the CNT content 
even when the CNT content is over 10 wt%. 



Verification 

To verify the effectiveness of the above multi- scale nu- 
merical simulations, the following theoretical prediction 
and experimental measurements were carried out. 

Theoretical prediction 

The following assumptions are made to derive conven- 
tional micromechanics models for the coefficient of ther- 
mal expansion (CTE). Note that the CTE, which is 
generally understood as a constant and temperature-inde- 
pendent, is different from the thermal expansion rate used 
here. Following the terminology of conventional microme- 
chanics models, we still use CTE in this section. The two- 
phase composite consisting of matrix and short fiber is of 
perfect interfaces at phase boundaries. Therefore, it is im- 
possible for the two components, i.e., the matrix and short 
fiber, to separate at their interfaces when the composite is 
loaded or heated. Additionally, only macro-composites are 
considered, namely the scale of the reinforcement is large 
compared to that of the atom size or grain size so that 
composite properties can be modeled by continuum 
methods. This assumption may be reasonable here since 
the present MWCNT is comparatively large in diameter. 
Finally, the composite properties are an appropriate aver- 
age of those of the components. 

The CTE of a composite with short-fiber orientation 
distribution function /((/?), which is independent of di- 
mension, can be given by [18] 



matrix. If n is the number of possible orientations, the 
CTE of the nanocomposites is 



VfE/af + V m E m a n 
VfEf + V m E m 



f 

Jo 



( cos 2 (p — v c sin 2 xp)f '((f) dip 



+ 



(l + v m )a m V m +(l + v f )a f V f 



f((f)dip. 
(1) 



For nanocomposites which contain a uni-directionally 
aligned reinforcement phase (e.g., MWCNT), f(cp) = 1, 
and therefore, the CTE of the nanocomposites is 



V/Ef a/ + V m E m a n 
VfEf + V m E m 



(2) 



If MWCNTs are randomly orientated, the orientation 
distribution function flip) = II n, where n represents the 
number of different orientations of the MWCNTs in the 



If VfEfCCf + V m E m a m 

2\ V f E f + V m E m [ ~ Vc) 

Hl + v m )a m V m +(l + v f )a f V f y 



(3) 



In the above equations, the nomenclatures for the 
parameters are as follows: 

a, CTE 

V, volume fraction 
E, Young's modulus 
v, Poisson's ratio 

and the subscripts are as follows: 

c, nanocomposite 
m, the matrix 

f, the reinforcement phase (MWCNT here) 

Note that Poisson's ratio of the nanocomposites, v c in 
Equation 3, was directly obtained from the rule of mix- 
ture and the data in Table 2. For 1-5 wt% addition of 
CNTs, v c ranges from 0.338 (1 wt%) to 0.333 (5 wt%). 

Experimental measurements 

In the present experiments, MWCNTs were made via 
chemical vapor deposition, with purity above 99.5% 
(Hodogaya Chemical Co., Ltd., Tokyo, Japan). The detailed 
data have been listed in Tables 1 and 2. An insulating 
bisphenol-F epoxy resin (JER806, Japan Epoxy Resins Co., 
Ltd., Tokyo, Japan) and an amine hardener (Tomaido 245- 
LP, Fuji Kasei Kogyo Co., Ltd., Osaka, Japan) were used as 
matrix. The MWCNT/epoxy nanocomposites were pre- 
pared by mixing the epoxy and the hardener using a 
planetary mixer (AR-100, THINKY Co., Ltd., Tokyo, 
Japan) at 2,000 rpm for 30 s. Then, the MWCNTs were 
added into the mixture and mixed again at 2,000 rpm for 
10 min. The final mixture was poured into a silicon mold 
and cured in a vacuum oven at 80°C for 2 h. This nano- 
composite fabrication method was the same with that in 
the authors' previous experimental work [19-21], in which 
very good dispersion states of the MWCNTs under 3 and 
5 wt% loading were identified (see image from scanning 
electron microscope observation in Figure 8 for the frac- 
tured surface of a 3 wt% sample). 

The thermal expansion properties of the MWCNT/ 
epoxy nanocomposites were measured using a TMA 
equipment (TMA-50, Shimadzu Co., Kyoto, Japan). The 
TMA measurement methodology is described as follows: 
a rectangular sample (3 cm wide, 3 cm long) was cut 
from the nanocomposites at a point 3 cm from the 
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parallel portion of the tensile test specimen (according 
to JIS K 7197 [22]). Specimens were heated from 30°C to 
120°C at a scanning rate of 5°C/min in air for continu- 
ous measurements. The thermal expansion properties of 
pure epoxy were similarly measured for the same speci- 
men size and test conditions. Note that the highest test 
temperature, i.e., 120°C, is close to the glass transition 
point of bisphenol-F epoxy resin, which usually ranges 
from 120°C to 130°C, depending on fabrication condi- 
tions. In our tests, it was found that even at 120°C, the 
obtained thermal expansion rates were still normal and 
a molten or rubber-like state in epoxy was not identified. 

Comparison 

Figure 9 shows the comparison between the thermal ex- 
pansion properties of the MWCNT/epoxy nanocompo- 
sites as determined by multi-scale numerical simulations, 
theoretical analysis, and experimental measurement. In 
Figure 9a, for uni-directional models, the comparison be- 
tween the thermal expansion properties by multi-scale nu- 
merical simulation and theoretical prediction was given, in 
which the relative difference is lower than 15% for the 
results. In Figure 9b,c, for multi-directional models, the 
comparisons of experimental, simulated, and theoretical 
results were shown for different CNT contents (i.e., 1 and 
3 wt%). It can be found that the multi-scale numerical 
simulation results possess a similar trend to the theoretical 
prediction and experimental measurement as temperature 
increases. It should be noted that the relative difference is 
also lower than 15% for all three results. This implies that 
the present multi-scale numerical simulation is effective in 
predicting the thermal expansion properties of CNT- 
based nanocomposites under the condition that the CNT 
is of a comparatively large size and a good dispersion state 
in matrix. Figure 10 shows the influence of CNT loading 
on the thermal expansion rates of the MWCNT/epoxy 
nanocomposites at high temperature (120°C), which was 
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Figure 9 Comparison of experimental, numerical, and theoretical 
results, (a) Simulated and theoretical results (uni-directional 
CNT/epoxy nanocomposite), (b) experimental, simulated, and 
theoretical results for 1 wt% (multi-directional CNT/epoxy 
nanocomposite), (c) experimental, simulated, and theoretical results 
for 3 wt% (multi-directional CNT/epoxy nanocomposite). 



evaluated by experimental, simulated, and theoretical 
approaches. From this figure, it can be found that the ther- 
mal expansion rate obtained by experiments decreases 
about 25% at 1 wt% and 35% at 3 wt%. Moreover, a similar 
trend is observed at a broad temperature range from 30°C 
to 120°C, in which the thermal expansion rate decreases 
with CNT loading for each case, and the present 
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Figure 10 Relationship between CNT content and thermal 
expansion rate of CNT/epoxy nanocomposite at 120°C. 



numerical simulation and theoretical analysis can effect- 
ively predict the experimental measurements. This indi- 
cates that the addition of CNT leads to a considerable 
reduction in the thermal expansion rate of the MWCNT/ 
epoxy nanocomposites, where CNTs with comparatively 
large size are well dispersed in the matrix in the present 
study. This characteristic leads to some special potential 
applications, such as good dispersion of CNTs into the 
matrix of carbon fiber-reinforced plastic to reduce residual 
stresses induced in the fabrication process. However, in 
many practical experiments, both distribution and disper- 
sion of the CNTs may be nonuniform because of the differ- 
ent properties of CNTs and fabrication methods; practical 
agglomeration of CNTs in the matrix may weaken this 
positive effect, i.e., reduction of the thermal expansion rate 
of the matrix. 

Conclusions 

In this work, the thermal expansion properties of CNT/ 
epoxy nanocomposites with CNT content ranging from 1 
to 15 wt% were investigated using a multi-scale numerical 
technique in which the effects of two parameters, 
temperature and CNT content, were investigated exten- 
sively. For all CNT contents, the obtained results clearly 
revealed that within a wide low-temperature range (30°C ~ 
62°C), the nanocomposites undergo thermal contraction, 
and thermal expansion appears in a high-temperature range 
(62°C ~ 120°C). It was found that at any CNT content, the 
thermal expansion properties vary with the temperature. As 
temperature increases, the thermal expansion rate increases 
linearly. However, at a specified temperature, the absolute 
value of the thermal expansion rate decreases nonlinearly 
as the CNT content increases. Moreover, the results pro- 
vided by the present multi-scale numerical model are veri- 
fied with those obtained from a micromechanics-based 
theoretical model and from experimental measurement. 
Therefore, this multi-scale numerical approach is effective 



to evaluate the thermal expansion properties of any type of 
CNT/polymer nanocomposites. 
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